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ABSTRACT 

Observations of distant sources of high-energy (HE) 7-rays are affected by attenuation resulting from the interaction 
of the 7-rays with low energy photons from the diffuse meta-galactic radiation fields at ultraviolet (UV) to infrared 
(IR) wavelengths (Extragalactic Background Light; EBL). Recently, a large data-set of HE observations from the 1st 
year survey of the Large Area Telescope (LAT) instrument on-board of the Fermi satellite became available, covering 
an energy range from 100 MeV up to 100 GeV. In this paper, the potential of such large HE data-sets to probe the 
density of the EBL - especially in the UV to optical - is explored. The data from the catalog is investigated for an 
attenuation signature in the energy range 10-100 GeV and the results are compared with the predictions from EBL 
model calculations. No clear signature is found. The statistics are still limited by (1) the sensitivity of Fermi/LAT to 
detect sources above 10 GeV, (2) the number of firmly identified sources with known redshift, both which will improve 
over the coming years. 
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1. Introduction 

Measuring the energy spectrum of distant sources of high- 
energy 7-rays (100 MeV - 100 GeV; HE) provides a unique 
method to probe the densities of extragalactic diffuse pho- 
ton fields at ultraviolet to infrared wavelengths (extragalac- 
tic background light; EBL). The HE 7-rays interact with 
the low energy photons of the EBL via the pair-production 
process leavin g an energy dependent attenuation sign ature 
in the spectra ( |Nikishov|1962| [Gould fc Schreder|1967| ). The 
strength of this attenuation and its spectral behavior de- 
pends on the properties of the EBL (i.e. its spectral number 
density) . 

Very-high energy (VHE; E>100GeV) 7-ray observa- 
tions of distant sources performed with ground-based in- 
strumen ts have successfully been used to limit the EBL den- 
sity (e g.lAharonian et al.|2006|[Mazin fc Raue|2007||Albert| 
|et al.||2008| ). Unfortunately, typical ground-based VHE in- 
struments have limited observation time (~ lOOOh/year) 
and field of view (2.5 — 5° diameter) and therefore the num- 
ber of known extragalactic VHE sources is small (O(30)). 
EBL studies using VHE data, therefore, often rely on the 
observations of a few sources (sometimes even a single 
source). 

Different possibilities ari se from observations with the 
Large Area Telescope (LAT; |Atwood et al.||2009[ ) on board 
the Fermi satellite. The detector has a large field of view 
(2.4 sr) and is operated mostly in survey mode covering the 
full sky every 3h. Fermi/LAT observations cover the en- 
ergy range from 100 MeV to 100 GeV and thereby enable 
to probe the EBL at ultraviolet to optical wavelengths. 
For the redshift range investigated in this paper (up to 
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z = 2 — 3), the energy spectrum below 10 GeV is essen- 
tially not affected by EBL attenuation. Strong contribu- 
tions from early stars to the EBL density could, in principal, 
lead to some attenuation in this energy range (Kashlinsky 
2005b ) but the claimed ex perimental evidence for such a 
high contribution (see e.g. |Dwek fc Arendtp998 Gorjian 



et al.||2000| [ Cambres yet aL||2001| JMatsumoto et al. 2005 
and Kashlinsky 2005 a] for a review) has been challenged on 



theoretical (e.g.|Madau fc Silk|2005) and experimental (e.g. 
Dwe k et al. |2005| |Aharonian et al . 2006; Thomps on et al. 
2007) grounds. Therefore, it is possible to (i) sample parts 
of the intrinsic spectrum emitted at the source and (ii) to 
observe sources up to higher redshift compared to VHE ob- 
servations^] With the recent rele ase of the Fermi/LAT 1 
year point source catalog (1FGL; Abdo & for the Fermi- 
LAT Collaboration|2010#] CAT1 in the fo llowing) and the 
T AGN catalo g (1LAC; [Abdo fc for the 
Fermi-LAT Collaboration! |2010c| CAT2 in the following) 



a large high-energy 7-ray data-set became available. The 
catalog contains ~700 7-ray sources associated with extra- 
galactic objects and enables, for the first time, to perform a 
population study for the global EBL attenuation effect on 
a sufficiently large data-set (e.g. Chen et "aT||2004 ) . 

In this paper the 1FGL data-set will be examined for 
a signature from EBL attenuation. The methods used for 
the study are discussed in Sec. [2| results are presented in 
Sec. [3] and Sec. [4] gives a summary and conclusions. For 
the calculations in the paper the following values for the 



1 Source with a redshift of z > 0.5 are expected to be severely 
attenuated due to the interaction with the EBL at energies above 
100 GeV. 

2 http:/ /fermi. gsfc.nasa.gov/ssc/data/ access/lat/lyr_catalog/ 
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Fig. 1. Energy dependent attenuation resulting from the 
interaction of HE/VHE 7-rays with the low energy photons 
from the EBL for differen t redshifts z (EBL model: KN low 
|2QQ2 ). The grey shaded area gives the 
highest energy bin in the 1FGL catalog 



from Kneiske et al 
energy range of the 
which will be investigate for attenuation features. 



cosmological parameters are adopted: h — 
tt M = 0.3. 
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2. Methods and data sample 

The 1FGL catalog lists for each source the integral flux for 5 
bins in energy (0.1-0.3 GeV, 0.3-1 GeV, 1-3 GeV, 3-10 GeV, 
and 10-100 GeV) and the results of the fit of a power law 
to the energy spectrum over the full energy range from 0.1 
to 100 GeV. Fig. [I] dis plays the EB L atte nuation for the 
low EBL model from |Kneiske et aL (2002) for the energy 
range 1 GeV to 200 GeV for several redshifts. It can be seen 
that significant EBL attenuation is only expected for the 
flux in highest energy bin. The sources in the catalog will 
therefore be examined for an attenuation signature in the 
highest energy bin (i.e. a deviation from the behavior at 
lower energies). The main caveat for this type of approach 
is that intrinsic curvature towards the highest energies is 
expected for many emission models and has been detected 
in several sources for energies below 10 GeV. Therefore, no 
(strong) conclusion on the EBL attenuation can be drawn 
from the fact that the flux in the last bin is attenuated. On 
the other hand, if no attenuation is detected, this can be 
used to set limits on the attenuation strength and therefore 
the EBL density. To search for a possible attenuation sig- 
nature in the highest energy bin two different methods will 
be employed: 

1. Power law fit. For this method, the power law fit to the 
energy spectrum is assumed to be a good representation 
of the overall spectrum emitted at the source. Due to 
the shape of the energy spectra^] and the sensitivity of 
the detector (best sensitivity in the 0.1-0.3 GeV energy 
range) the spectral fit result for the full energy range 
are dominated by the photons with energies < 10 GeV. 
The integral flux in the energy range 10 to 100 GeV as 
calculated from the power law fit Fi nt pL can therefore 
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Table 1. Number of sources in the 1FGL catalog for spe- 
cific criteria: i^io-iooGeV - significant flux in the highest 
energy bin; z - known redshift; NC - no significant cur- 
vature detected in the spectrum; C - significant curvature 
detected in the spectrum; FSQR - source is of type FSRQ. 



be considered a good representation of the flux for the 
un-attenuated spectrum. The ratio of the measured flux 
in the last bin i^io-iooGeV and Fi nt pL is used as an 
indicator to quantify the strength of the attenuation. 
The main caveat for this approach is that it relies on 
the power law fit to be a good representation of the 
overall spectrum. 
2. Flux stacking. The flux in the highest energy bin 
^10-100 GeV is compared to the integrated flux at en- 
ergies <10GeV Fo. 1-10 GeV by summing up the flux 
in the bins 0.1-0.3 GeV, 0.3-1 GeV, 1-3 GeV, and 3- 
10 GeV. The ratio between the two integral fluxes 
F10-100 GeV /Fo. 1-10 GeV is used as indicator for an at- 
tenuation signature. This ratio depends on the spectral 
shape of the source and therefore, to be able to compare 
and combine the results for different sources, a sample 
of intrinsically similar sources is needed. While the BL 
Lac type sources show a wider spread of spectral shapes 
(spectral indices), the spectral shapes of the flat spec- 
trum radio quasars (FSRQs) in the Fermi/LAT range 
are remarkably similar, displaying a narrow distribu- 
tion of spectral indices with a mean of T ~ 2.5 (see 
CAT2 Fig. 12). Therefore, only spectra from FSRQ type 
sources will be used for this type of analysis. 

EBL attenuation is a redshift dependent process. A sig- 
nature should therefore emerge when analyzing the evo- 
lut ion of the a ttenuation effect with redshift (for caveats 
Reimer 2007). The 1FGL catalog does not contain red- 



3 Most AGN energy spectra in the 1FGL are well described 
either by power laws with T > 2 or log parabola spectra. 



shifts, but these can be found (for the associated sources) 
in the 1LAC catalog. For the analysis the two catalogs are 
merged based on the catalog names using only the high 
quality data sample from the a 1LAC catalog (clean sam- 
ple, see CAT2). Of the 1451 sources in the combined cat- 
alog 287 (~ 20%) have a significant detection in the last 
energy bin, 390 (~ 27%) have a redshift associated, and 72 
(~ 5%) have both. 57 of these source do not show significant 
curvature (curvature index < 11.34, see CAT1), leaving 15 
sources with significant curvature in the spectrum. For the 
2nd method, only FSRQ type sources will be used, which 
do account for 12 sources without and 6 sources with sig- 
nificant curvature. All numbers are summarized in Tab. Q] 
The number of sources used in this study is only a small 
fraction of the total number of sources in the 1FGL. How 
this impacts the results and what to expect in the future 
will be discussed in Sec. [3] and Sec. [4] respectively. 
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Fig. 2. Ratio of the integral flux in the highest energy bin 
^io- ioo GeV and the flux expected from the power law fit 
FintPL versus redshift for AGN sources from the 1FGL cat- 
alog (method 1). For comparison the expectation for differ- 
ent EBL models is plotted. Error bars are la and include 
the uncertainty due the uncertainty in the power law fit. 



3. Results 

Power-law fit. Figure [2] shows the ratio Fio-ioo GeV /FintPL 
versus redshift for tne source samples described in the 
previous section. For comparison, the expectation for the 
EBL attenuation for three different EBL models is also 



shown ( |Franceschini et aL 2008 = Fran08 ; KN low/high = 
low/high model from Kneiske et al.||2002 ). The models are 
chosen to cover some range in attenuation with the Fran08 
model being at the low end, describing the EBL density 
derived from galaxy counts, and the KN high model re- 
sulting in the strongest attenuation but overproducing the 
results from galaxy counts. The integrated attenuation for 
the EBL models is calculated under the assumption that 
the intrinsic spectrum follows a power law with spectral in- 
dex T = 2.5 j For redshift s < 0.5, where the majority of the 
sources in the sample is located (63%), the expected EBL 
attenuation signature is less then 10%. For redshift z = 2.5 
up to 50% attenuation is expected for i^io-iooGeV- For the 
sources without intrinsic curvature the errors on the flux 
ratio are of the order of 0.2 to 0.5. The errors include the 
uncertainty in the parameters of the power law fit. At lower 
redshifts the flux ratios scatter around 1 as expected. At 
higher redshifts (z > 1.25) only very few sources are avail- 
able (4) and no clear attenuation signature is visible (i.e. 
sources following the trend of the theoretical predictions). 
For completeness, the results for the sources with curved 
spectra are also shown. The curvature affects the power 
law fit and therefore the absolute value of the attenuation 
has to be taken with caution since it is likely underesti- 
mated. As expected, the sources with significant curvature 
show deviation from the power law fit for i^io-iooGeV- For 
redshift > 0.5 all curved sources show a concave behavior 
(attenuated last bin) and they all lie below the expectation 
for the EBL model with the highest attenuation used in 
this study 



FSRQ curvature 
Mean z lo Jz high 
EBL model (Fran08) 
EBL model {KN high) 



4 The change in integrated attenuation when using Y 
" = 3 is less then 10%. 
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Fig. 3. Ratio of the fluxes in the energy range 10-10 GeV 
and 0.1 to 10 GeV i^io-iooGev/^o.i-ioGW versus redshift 
for FSRQs from the 1FGL catalog (upper panel: without 
curvature; lower panel: with curvature). In blue the mean 
values for two bins in redshifts are shown (z < 1.25 and 
z > 1.25; filled symbols: flux stacking; open symbols: av- 
erage; see text for details). The colored bands denote the 
expectation for two different EBL models normalized to the 
low redshift mean value. In the upper panel one point at 
Z = 0.22 With Fio-iOOGev/^0.1-10GeV = 0.014 ± 0.007 is 
not shown to improve the plots readability. 



Flux stacking. The Fio-iooGev/^o.i-iOGeV values from the 
flux stacking for FSRQs without and with significant cur- 
vature are shown in Fig. [3] upper and lower panel respec- 
tively. Also shown are the mean values for two bins in red- 
shift (z < 1.25 and z > 1.25). The mean is calculated us- 
ing two different methods: (1) The Fio-iooGev/^b.i-iOGeV 
values for different sources in one redshift bin are aver- 
aged, weighted with 1/a 2 (open blue squares). (2) For all 
sources in one redshift bin, first, the fluxes are averaged in 
the different energy bins (again weighted with 1/cr 2 ) and 
then summed to get the total fluxes for the ranges 0.1-10 
and 10-100 GeV (filled blue squares). The second method 
has the advantaged that sources, which do not have a sig- 
nificant detection in one of the energy bins, do not dis- 
tort the summed flux. The mean values are plotted at the 
mean redshift of the sources in the bin. For comparison, 
the expectations for different EBL models are shown. Since 
^10-100 Gev/Fo. 1-10 GeV is a relative value the EBL model 
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predictions are normalized to the mean value in the low red- 
shift bin (at the mean redshift). Shown is the outer value of 
the one sigma error band for the two different mean values. 
As noted earlier, the total number of sources in the analysis 
is low, in particular for the case of the curved spectrum (6 
sources in total). For the individual sources the spread is 
rather large and no clear trend is visible. 

When considering the mean values for the two bins in 
redshift it can be noted that, for the source with no curva- 
ture (upper panel), the mean in the high redshift bin lies 
above the value for the low redshift bin. While this is not 
significant it might indicate a change in the intrinsic spec- 
tral properties with redshift. The mean value for the high 
redshift bin lies also ~ 2<j above the prediction for the high 
EBL model (KN high). Higher statistics and a closer look 
at the individual sources, which is beyond the scope of this 
paper, will be needed to investigate this further. 



4. Summary & Conclusion 

In this paper the HE data from the 1FGL catalog is in- 
vestigated for an attenuation signature resulting from the 
interaction of HE 7-rays and low energy photons from the 
EBL. Such an EBL attenuation signature is expected to 
show up in the energy range 10-100 GeV and should exhibit 
a characteristic redshift dependence. The catalog provides 
a homogenous sample of HE data which, for the first time, 
enables a population study on the EBL attenuation effect 
on such a sample. Two different methods have been applied 
to search for such a signature but no significant signal has 
been detected. The expected signal is rather small and the 
current statistic in the catalog is not sufficient for a strong 
detection with the methods applied. 

The 1FGL contains results from 11 month of data tak- 
ing. Fermi is expected to operated for five to ten years which 
will significantly increase the statistics. The error on the 
power law fit and on the flux in the last bin Fio-iooGeV 
will improve by a factor y/b ~ 2 to y/lO ~ 3. The in- 
creased statistic will also lead to the significant detec- 
tion in the highest energy bin for many of those sources 
for which at the moment only an upper limit is avail- 
able. The number of sources having a significant detec- 
tion in the highest energy bin after 5/10 years of obser- 
vations can be estimated using the results from the power 
law fit given in the catalog: adopting a minimum flux of 
^10-100 GeV — 1-5 • 10 _10 cm _2 s _1 for a significant detec- 
tion in the 11 month data-set0 it can be estimated that 
after 5 (10) years of observations ~50% (~65%) of the 
AGN sources in the 1FGL will have a significant detec- 
tion at 10-100 GeV compared to 20% after 11 months. The 
increased observation time will also lead to the discov- 
ery of new sources not yet included in the 1FGL. Of the 
~700 sources in the 1FGL associated with extragalactic 
sources only ~55% have a known redshift. Dedicated multi- 
wavelength observation campaigns are being performed to 
improve the data available for the associated sources (e.g. 
redshift). In addition, the understanding of the systematics 
at the high energy end of the Fermi/LAT energy regime, 
which, for the moment, has possibly still large uncertain- 
ties, will greatly improve. 



Additional information can be derived by combin- 
ing the Fermi/LAT HE data with data from lower and 
higher energies. Simultaneous observations at low energies 
(radio/optical/x-lrays) enable to constrain the modeling of 
the sources and thereby enables to de rive constraints on 
the high energy part of the spectrum ( Mankuzhiyil et al. 
|2010[ ). Combining the HE data with VHE observations from 
ground-based instruments gives a han dle on the spectrum 



with and without EBL attenuation (Abdo et al. 
although this type of study is limite< 



2009 ) 

by the number of 



sources detected with both instruments and the difference 
in sensitivity. This situation will improve with the upcom- 
ing generation of ground-based VHE instruments (MAGIC 
II, H.E.S.S. II) which focus on the energy regime between 
10 and 100 GeV. 

Note: After submission of this article (May 3, 2010) an ar- 
ticle b y the Fermi collaboration appeared on the preprint 
server ( I Abdo fc for the Fermi-LAT Collaborationist) 10b ) 
which discusses, among other methods, a similar investi- 
gation to search for an EBL attenuation signature in the 
Fermi catalog data, which gives comparable results. 
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